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Complex photonic band structures (CPBS) of transmission metallic gratings with rectangular 
f*^*; | slits are shown to exhibit strong discontinuities that are not evidenced in the usual energetic band 

structures. These discontinuities are located on Wood's anomalies and reveal unambiguously two 
different types of resonances, which are identified as horizontal and vertical surface-plasmon reso- 
nances. Spectral position and width of peaks in the transmission spectrum can be directly extracted 
i— ( ■ from CPBS for both kinds of resonances. 
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■ Recently, metallic films with sub-wavelength apertures became the subject of increasing interest. Experiments show 
that 2D arrays of sub- wavelength holes in metallic films can lead to extraordinary transmission of light . These 
properties have been attributed to the excitation of coupled surface-plasmons on the upper and lower surfaces of the 
grating. The ability of these structures to control light has been shown [||J7), and applications have already been 
proposed in order to exploit these properties in different fields, such as electromagnetic filters or photolithography. 
£L| Electromagnetic calculations have been carried out by Porto et al. || in order to study the mechanisms that 
enhance the transmission of light through metallic gratings with very narrow slits. They distinguished two different 
CO ■ mechanisms, that is the excitation of coupled surface plasmon polaritons on both surfaces of the metallic grating, and 
the coupling of incident plane waves with waveguide resonances located in the slits. They correspond respectively to 
surface-plasmon bands and flat bands in the energetic band structure. Similar resonances located in grooves have also 
been observed experimentally and numerically in the reflected light of metallic gratings |p|-|IT| . 

Complex photonic band structures (CPBS) were previously calculated by Kuzmiak et al. for 2D periodic systems 
with metallic components p2j| , in order to obtain the attenuation and the lifetime of each mode. CPBS allowed a 
splitting of the lifetime of degenerate modes at Brillouin-zone boundaries to be observed. 

In this Letter, we will show that, in the case of rectangular metallic gratings, complex dispersion curves demonstrate 
the existence of a new kind of discontinuities. Much stronger than for 2D periodic systems, they are located on Wood's 
anomalies. They correspond to the transition between two different types of resonances, identified as horizontal and 
vertical surface-plasmon resonances. The first one is a periodic structure resonance, the second one is a Fabry- 
, Perot like resonance. The discontinuity can be up to four orders of magnitude. Moreover, we will see on zero-order 
transmission spectra that the calculated complex frequencies are in excellent accordance not only with the spectral 
| position, but also with the width of resonances' peaks. 

Fig. [j] shows the structure studied in this Letter, which is similar to the one of Porto et al. ||. In the following, 
the period of the grating d = 3.5 /im and the width of the slits a = 0.5 /xm will be kept constant, whereas the height 
h of the grating will vary. It is a symmetric structure, surrounded with air. The metal is assumed to be gold, whose 
dielectric function is taken from the reference used by Porto Ejj ]. 

Our calculations were carried out for the TM polarization (the incident magnetic field is parallel to the direction y 
of the grating, see Fig. [I]). They are based on the exact modal expansion in region II |l4j and S*-matrices formalism. 
In the upper and lower regions (I and III), the electromagnetic field is expressed as superposition of plane waves 
e Mfc<"'i±t< T,) z)-jwt] (Rayleigh's expansion). Their pure real wave vector component fci™\ which represents the in-plane 
momentum of photons is given by Eq. (TH): 



X, 

* : 4 n) = ^ + (i) 

The two other terms and u> are complex quantities, according to the relation ^fci"^ + ^fci™^ = (lu/c) 2 . 

For the calculation of CPBS, we use a S'-matrix formalism, in which the amplitude of the reflected and transmitted 
diffracted waves are given by a matrix S applied to a vector containing the amplitude of the upper and lower incident 
waves. An electromagnetic mode is defined by the condition of the existence of a wave in the structure without any 
incident wave, which is equivalent to say that the inverse (5* -1 ) of the ^-matrix has an eigenvalue that is null, with 
an eigenvector (i? n ,T n ), corresponding to an electromagnetic field in regions I and III: 
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with 5R(i/ci n) ) < 0. 

For a given k x — k x °\ the electromagnetic modes are defined by complex frequencies loq. The inverse of the 
imaginary part of the frequency provides the amplitude lifetime 2t of the resonance fl5|,|l2| , which can be associated 
to the dimensionless quality factor Q — 5R(wo) t. 

The complex dispersion curves are lines in the three-dimensional space defined by the complex frequency loq of 
resonance modes and the parallel momentum of photons k x . For a more comprehensive representation we are plotting 
two projections of ujq, the energy E — % Sft(wo) and the dimensionless quality factor Q = 5ft(wo)/(~2 3(wo))- 

Fig. shows the complex dispersion curves calculated for a height h = 1.4 /im. The lines on the energetic band 
structure (Fig. ^(a)) represents the frontier of the Wood-Rayleigh's anomalies. They correspond to the emergence (or 
the vanishing) of a diffracted wave in the region I or III, which wave is then at the grazing angle. These anomalies 
are observed on transmission spectra when the n th parallel wave vector kx (given by Eq. (Q)) is equal to the real 
incident wave vector kg. 

The energetic band structure in Fig. g(a) shows two different types of modes, whose dispersion relations have already 
been studied by Porto. The bands located close to Wood's anomalies, with slightly smaller energies, characterize 
horizontal surface-plasmon resonances. They are excited by the first non-propagating (or evanescent) diffracted wave, 
and their quality factors are greater than 10 and up to more than 10 4 (see Fig. §(b)). On the other hand, flat bands 
represent electromagnetic modes which are almost independent of the incident angle 9, and with lower quality factors 
(below 20). They were previously reported as waveguide resonances We will justify further in this Letter why 
these resonances are more precisely vertical surface-plasmon resonances. 

The most remarkable result is given by the quality factor of the upper energetic band, which is broken into three 
parts. Considering the point A in Fig. ^(a), we sec that the energetic band is continuous whereas the quality 
factor's curve is discontinuous at the corresponding points A' and A" in Fig. |](b). On the left side of point A, the 
flat band corresponds to a vertical resonance with very low quality factors (of the order of 3), whereas the right side 
corresponds to an horizontal surface-plasmon resonance with high quality factors (up to more than 10 4 ). The lifetimes 
of the resonances are r = 2.9 fs at A' and r = 33 ps at A" . 

In order to understand physically the behavior of the electromagnetic waves, it must be pointed out that the 
discontinuous transitions between the two types of resonance are located exactly on Wood's anomalies (point A 
and its symmetrical point B), so that they are linked to the emergence or the vanishing of a diffracted wave. This 
phenomenon can be explained in considering the energetic band structure: below Wood's anomaly, the first non- 
propagating (evanescent in z direction) diffracted wave (n = +1 at the point A and n = -2 at the point B) is 
responsible for the excitation of the horizontal surface-plasmon resonance. When the resonance crosses Wood's 
anomaly, the evanescent wave becomes a z-direction propagating diffracted wave (at grazing angle) which can not 
excite an horizontal resonance anymore, but is coupled to vertical surface-plasmons located on the vertical walls of 
slits. Then, we obtain a flat band which corresponds to a "weaker" resonance in the Fabry-Perot like cavity delimited 
by the top and bottom surfaces of the grating. This vertical resonance is coupled with the last propagating diffracted 
wave, thus it ends at the crossing with Wood's anomaly, where this wave vanishes. 

The complex band structure in Fig. || has also a continuous transition between the horizontal and the vertical 
resonances in the lower band. The flat band at normal incidence corresponds to vertical resonances, and becomes a 
surface-plasmon band for greater incident angles. Its quality factor increases from 20 to about 300 during this smooth 
transition. The last flat band at about 0.35 eV corresponds to a vertical resonance for the whole range of k x in Fig. 

By varying the width a of the slits, we have observed that the complex mode of the vertical resonance is nearly 
unchanged. In the calculation, it can be seen from the modal decomposition of the electromagnetic field, that this 
resonance corresponds mainly to the excitation of the first eigenmode of a 2 infinitely extended region II. The wave is 
evanescent in the x direction and propagates at the vertical interfaces air/metal (z direction). This TM mode exists 
for any width of slit, and becomes equal to the wave vector k ps of the surface-plasmon of a semi-infinite air-metal 
interface in the limit of wide slits, so that it is more correct to call it a vertical surface-plasmon mode instead of a 
waveguide mode. The transitions between the regions II/I and II/III act as mirrors for the vertical surface-plasmon 
mode. Their reflection coefficients are low, so that the vertical resonance is a Fabry-Perot like resonance with low 
quality factors. 

Complex dispersion curves shown above can be seen as convenient compact representations of the resonances, which 
can be associated to transmission spectra. For that purpose, we make the hypothesis that the complex frequency lOq 
of a resonance i is a single pole of the amplitude of the transmitted waves. Consequently, the zero-order transmitted 
energy can be approximated by a Lorentzian function near the resonance i. Assuming u> is the (real valued) frequency 
of the incident wave, we obtain: 
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To = ■ (3) 

The complex band structure of the grating is shown in Fig. ^ for a height h — 3 jjm. Fig. |^ shows the zero-order 
transmission spectra calculated for the same height and two different angles of incidence. They correspond to 9 = 0° 
(a) and 9 = 12° (b) cuttings in Fig. ||(a). The Lorentzian curves corresponding to Eq. (||) are added to the spectra. 
The constant a, is fitted to obtain the same maximum value as the spectra. At normal incidence (Fig. |](a)), the 
broad resonances at 0.38 eV and 0.54 eV (not represented in Fig. ||) are not separated enough to allow the comparison 
to their respective Lorentzian functions to be made. In such a case, the amplitude of the transmitted wave is the sum 
of two complex Lorentzian functions, and phase considerations make it more difficult to fit to the spectra. 

For the other peaks, Fig. ^ shows an excellent agreement of the Lorentzian curves with the spectra. Hence, the 
single-pole approximation is justified and Eq. (^) gives the width at half maximum: 

Aw = -23K). (4) 

It justifies also the above definition of the quality factor of the resonance, which corresponds to the more conventional 
Q = E/AE. 

These transmission spectra allow some important characteristics of the resonances to be pointed out. At first, it 
must be noted that the flat band at about 0.35 eV has no solution at normal incidence, and the peak 7 disappears 
at 9 = 0°. An energetic shift can be seen for the horizontal surface-plasmon resonances between 9 — 0° and 9 = 12°, 
whereas the vertical resonance at about 0.17 eV does not move. Another point which comes out of the calculations 
is that vertical resonances in contiguous slits are nearly uncoupled. As a consequence, a periodic structure is not 
necessary for them to exist: the vertical surface-plasmon resonance could be excited in an unique slit, and fully coupled 
with a few ones. Hence, an important result is that the vertical resonance can be efficiently excited by a strongly 
focalized beam. 

In conclusion, we have carried out the calculation of the complex dispersion curves of rectangular metallic gratings 
in TM polarization. CPBS are directly linked to experimental data, like transmission spectra. The quality factor and 
the lifetime of resonances, that they provide, are crucial characteristics in many applications already proposed for the 
extraordinary effects appearing in such transmission metallic gratings. CPBS show pronounced discontinuities located 
on Wood's anomalies, which can not be evidenced in the usual energetic dispersion curves. They allow to distinguish 
unambiguously the two different natures of resonances: horizontal and vertical surface-plasmon resonances. Similar 
results can be obtained for reflection gratings. This work can also be extended in order to analyze resonances in 
more complex geometries, as Ebbesen's structures |lj for example, for which a cutoff diameter exists for vertical 
modes. The characteristics of the resonances studied in this Letter open up new possibilities in photonic materials. 
Moreover, these structures allow huge concentrations of light in very small volumes, which could have applications in 
fundamental optics Jig] and in optoelectronics. 
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FIG. 1. Description of the metallic grating. Horizontal and vertical surface-plasmons are symbolically represented. 



FIG. 2. Energetic (a) and quality-factor Q (b) representation of the first four photonic bands of a rectangular metallic 
grating (d = 3.5 fim, a — 0.5 (im) with a height h = 1.4 /jm. The lower band is only computed for energies higher than 0.125 
eV (A < 10 pm). Wood's anomalies (lines) are added to the energetic band structure (a). 

FIG. 3. Energetic (a) and quality-factor Q (b) representation of the first five photonic bands of a rectangular metallic 
grating (d — 3.5 /jm, a = 0.5 (im) with a height h — 3.0 /jm. The lower band is only computed for energies higher than 0.125 
eV (A < 10 pm). Wood's anomalies (lines) are added to the energetic band structure (a). 



FIG. 4. Comparison between zero-order transmission spectra (solid lines) and Lorentzian curves (dotted lines) for a 
rectangular metallic grating (d = 3.5 /im, a = 0.5 /im, h = 3.0 /im). The angle of incidence is (a) 9 = 0° and (b) 6 — 12°. 
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